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DN interaction from meson-exchange and quark-gluon dynamics 
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Abstract. We investigate the DN interaction at low energies using a meson-exchange model supplemented 
with a short-distance contribution from one-gluon-exchange. The model is developed in close analogy to 
the meson- exchange KN interaction of the Jiilich group utilizing SU(4) symmetry constraints. The main 
ingredients of the interaction are provided by vector meson (p, uj) exchange and higher-order box diagrams 
involving D* N, DA, and D* A intermediate states. The short range part is assumed to receive additional 
contributions from genuine quark-gluon processes. The predicted cross sections for DN for excess energies 
up to 150 MeV are of the same order of magnitude as those for KN but with average values of around 
20 mb, roughly a factor two larger than for the latter system. It is found that the w-exchange plays a 
very important role. Its interference pattern with the p-exchange, which is basically fixed by the assumed 
SU(4) symmetry, clearly determines the qualitative features of the DN interaction - very similiar to what 
happens also for the KN system. 
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1 Introduction 

The study of the interactions of charmed hadrons with 
nucleons is of interest in several contexts. One example is 
in experiments of relativistic heavy ion collisions (RHIC). 
Since long time [1] the suppression of J/W production in 
RHIC is being considered as a possible signature for the 
formation of a quark-gluon plasma (QGP). The alleged 
suppression would occur because the deconfined quarks of 
the QGP screen the long range confining potential thus 
making impossible the formation of the mesonic bound 
states. However, collisions of the charmed mesons with 
hadrons in the medium can also lead to dissociation of 
these mesons, subverting therefore the screening scenario. 
Moreover, more recently it has been argued that heavy 
quarkonia could be re-formed via rescattering processes 
of open-charm hadrons in the late stages of RHIC which 
would then lead to an enhanced J/& production [2]. Thus, 
it seems clear that a good knowledge of the interaction of 
charmed mesons with ordinary hadrons like nucleons is 
a prerequisite for differentiating between these scenarios. 
For a recent review on these issues, see Ref. [3]. Another 
example where the interaction of charmed mesons with 
ordinary hadrons is of interest refers to studies of chiral 
symmetry restoration in a hot and/or dense medium. In 
this respect, the interaction of charmed D mesons - which 
are composed of one light and one heavy quark - with 
nucleons is of particular interest. The properties of the 



light quarks in a D meson are sensitive to temperature 
and density and, therefore, changes in the properties of 
the D mesons in medium can be expected. Consequently, 
one can also expect that their interactions with nucleons 
will change in the medium. The D meson and its lowest 
excitations are somewhat special in this respect because 
their spectroscopy is simpler than of ordinary mesons com- 
posed solely by u and d quarks. This is so because the 
charm quark c is much heavier than the light u and d 
quarks, and to a good approximation these mesons can be 
described as one-body bound states, a fact that simplifies 
tremendously their study. 

Before one can infer in a sensible way changes of the 
interaction in the medium, a reasonable understanding of 
the interaction in free space is required. However, here 
one has to cope with a major difficulty, namely the com- 
plete lack of experimental data at low energies for the 
free-space interaction. This situation is hopefully going to 
change soon with the operation of the FAIR facility at the 
GSI laboratory in Germany. There are proposals for ex- 
periments by the PANDA collaboration [3] at this facility 
to produce D mesons by annihilating antiprotons on the 
deuteron and, through the rescattering of the produced D 
and D mesons on the spectator nucleon [5], to determine 
DN as well as DN cross sections and possibly even phase 
shifts. Still, for the design of detectors and of efficient data 
acquisition systems, estimates for the magnitude of such 
cross sections are urgently required. Therefore, there is a 
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need for developing models of the interaction of charmed 
particles with ordinary hadrons, and since - as said - 
not much is known empirically, such models can only and 
should be constrained as much as possible by symmetry 
arguments, analogies with other similar processes, and the 
use of different degrees of freedom. 

The interaction of charmed mesons with ordinary had- 
rons composed of u and d quarks has been investigated 
using effective hadronic Lagrangians and quark models. 
Most of the studies have concentrated on the interaction 
of J/W and other heavy charmonia with ordinary hadrons, 
mainly due to the interest in the QGP suppression hypoth- 
esis alluded above - see Ref. [6] for a review on these inves- 
tigations. With respect to the interaction of the D meson 
with the nucleon, which is the subject of the present pa- 
per, not much is known. The work of Ref. 7 , using an 
effective SU(4) hadronic Lagrangian, to the best of our 
knowledge was the first one to provide estimates of cross 
sections for the DN system in Born approximation. In 
terms of quark degrees of freedom, the authors of Ref. [8] 
have made an estimate for the DN cross sections using 
quark rearrangement arguments and concluded that such 
cross sections should be equal to the corresponding KN 
cross sections, though no explicit model was employed. 
(See also the results presented in Ref. [9]). 

In the present paper we investigate the DN interaction 
within a meson- exchange model and a quark model uti- 
lizing one-gluon-exchange (OGE), in the spirit of a recent 
study of the KN system by us |T0] . {D is used here generi- 
cally for the D° and D~ isospin doublet which contains a c 
quark, and corresponds to K consisting of the K + and K° 
isospin doublet with an s quark.) To be more specific, the 
DN interaction we construct is an extension of the KN 
meson-exchange model of the Jiilich group [TT 1 fT2 l [T3] , gen- 
eralized by assuming as a working hypothesis SU(4) sym- 
metry constraints. Note that the KN model described in 
Refs. [TT|ri2] considered not only single boson exchanges 
(cr, p, to), but also contributions from higher-order dia- 
grams involving N, A, K and K* intermediate states. We 
focus on the DN system because it has the advantage 
that its dynamics should be governed predominantly by 
the same "long-range" physics as the KN interaction, i.e. 
by the exchange of ordinary (vector and possibly scalar) 
mesons. Thus, fairly reliable and, most importantly, es- 
sentially parameter-free predictions can be made once one 
accepts the constraints provided by SU(4) symmetry. The 
DN system is expected to exhibit a much richer structure 
[14] and thus may be more interesting [T5 l lT6 l fT7] . but it 
involves also much larger uncertainties. In this case, like 
in the analogous KN system, there are couplings to sev- 
eral other channels which are already open near the DN 
threshold (A c tt, 2J c tt) or open not far from the threshold 
{A c rj). It is obvious that the coupling to those channels 
must play a crucial role for the dynamics of the DN sys- 
tem - as it does in the corresponding KN system - and, 
thus, will have a strong impact on any quantitative results. 
But the transitions to those channels and the interactions 
in those channels involve charmed baryon resonances as 
well as the exchange of charmed mesons, for example the 
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Fig. 1. Meson-exchange contributions included in the DN in- 
teraction. 



Z?*(2010), whose coupling constants and associated vertex 
form factors, required in any meson-exchange model, are 
practically unknown and difficult to constrain. 

Our DN model is obtained by substituting the one- 
boson-exchange contributions, but also the box diagrams 
involving K*N, KA, and K* A intermediate states, of the 
original KN model of the Jiilich group by the correspond- 
ing contributions to the DN interaction under the con- 
straint of SU(4) symmetry. Regarding the short-ranged 
quark part, we use the dominat contributions of OGE ex- 
change, which are the Coulomb and spin-spin parts. These 
are of the same form as for the KN interaction, but the 
mass of the charm quark is much heavier than the one of 
the strange quark, and the size parameter of the meson 
wave function is also different. Therefore, the DN inter- 
action will be different from the KN interaction. We want 
to emphasize that we iterate the effective meson-exchange 
and OGE potentials in a Lippmann-Schwinger equation, 
contrary to the common practice of using only Born ap- 
proximation 7,8,9. We found that the effect of iteration 
on the predicted cross sections can be quite substantial, 
being of the order of 50 % in some cases. 

The plan of our paper is the following. In the next sec- 
tion we specify the three-meson vertices used in the paper. 
Specifically, we discuss the constraints of SU(4) symme- 
try to relate the required new couplings and present the 
numerical values of vertex parameters. The interaction 
Lagrangians, which are needed to complete the deriva- 
tion of the meson-baryon potential, are presented in Ap- 
pendix [X] In Section [3] we discuss the quark-gluon ex- 
change mechanism for the DN interaction. The detailed 
derivation of the equations shown in this Section are out- 
lined in Appedix 151 Our numerial results are presented in 
Section |4j The paper ends with a short summary. 



2 The DN interaction in the meson-exchange 
picture 

The meson-exchange model of the DN interaction is con- 
structed in close analogy to the corresponding KN poten- 
tials developed by the Jiilich group some time ago [H1IT2] . 
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In those models, derived within time-ordered perturba- 
tion theory, not only single-meson (and baryon) exchanges 
were taken into account, but also higher-order box dia- 
grams involving K*N, KA, and K * A intermediate states. 
Thus, we will consider the corresponding contributions to 
the DN interaction too, cf. Fig. [TJ The general scheme 
and also the explicit expressions for the various contribu- 
tions to the interaction potential are described in detail in 
Refs. [TTiri2] and, therefore, we do not reproduce them in 
the present paper. We only summarize the used interac- 
tion Lagrangians in Appendix A. Here we want to focus 
on the SU(4) structure which is used to extend and relate 
the DN interaction to the KN system. 

For the construction of the DN interaction we need 
three-meson vertices involving charmed mesons of the kind 
PPV and VVP (P — pseudoscalar meson, V = vec- 
tor meson). The general form of the SU(4) invariant La- 
grangian is 

CmMM = g{15}[-aTr([M{ 15 y,M{ 15 y]M( 15 y) 

+ (l-a)Tr({M {15} ,M {15} }M {15} )} 
+ 3{i5}{i5}{i}(l - a)Tr({M {ls} , M {15} }M {1} ) 
+ 5{i5}{i}{i5}(l-a)Tr({M {15} ,M {1} }M {15} ) 
+ 9{1} (1- a)Tr({M {1} ,M {1} }M {1} ) , (1) 

where a is the F/ (F + D) ratio and M{ 15 } (Mm) stands 
for the SU(4) meson-15-plet (-singlet) matrix. For pseudo- 
scalar (P) and vector (V) mesons M{ 15 } are 4x4 matri- 
ces which are given in Table [T] Note that the PPV ver- 
tices involve only P-type coupling (a = 1) if we require 
charge conjugation invariance while the VVP vertices in- 
volve only P>-type coupling (a = 0). Let us stress again 
that SU(4) should not be a flavor symmetry of QCD but 
rather a working hypothesis to get a handle on the vari- 
ous couplings and form factors employed in our model. It 
is, of course, strongly broken due to the use of the very 
different physical masses of the various mesons. 

Based on the assumed SU(4) symmetry all relevant 
three-meson coupling constants can be derived from the 
empirically known ixixp coupling. In the Jiilich model [12j 
the value g nnp = 6.0 is used. The coupling constants of 
the other vertices that follow from this value are listed in 
Table [2] 

As far as the coupling constants belonging to the JVJV 
and NA vertices are concerned we take precisely the same 
values as in Ref. [T2], which are based on those of the 
(full) Bonn NN potential, cf. Ref. [IS]. These coupling 
constants are listed in Tabled] too. The Jiilich KN poten- 
tial contains also vertex form factors P that are meant to 
take into account the extended hadron structure and are 
parametrized in the conventional monopole or dipole form 
[TT1IT2] . In the present study of the DN system the cut- 
off masses appearing in those form factors for the various 
thrcc-mcson and baryon-baryon-meson vertices are like- 
wise taken over from Ref. [12] . Specifically, we make the 
assumption that Fj)Q m {q^) ~ FKKm(Qm)- This prescrip- 
tion is motivated by the notion that those form factors 
parametrize predominantly the off-mass-shell behaviour 



of the exchanged particles - which are indeed the same 
in the KN and in the DN interaction. 

Let us make some more comments about the coupling 
constants at the three-meson vertices. SU(4) symmetry 
implies the following for the vector meson coupling con- 
stants relevant for our study: 



VSgKKp = V3 - gTrrrp, 



9KKp, 9ddu 




Q'K'Kp 



9dd p = 9kk p - 2 



(2) 



(3) 



Assuming ideal mixing of the 0*15, us and u>\ one obtains 
for the coupling constants of the physical lj and <fi 



9ddu — \j ~^9dDu)! 



9dd4> 



9dduh 



9ddu 



9ddui 5 ■ ( 4 ) 



The same relation holds also for the K meson. In case of 
the K meson the coupling constant gxKu is given by that 
of 9kku> s alone, since there is no singlet coupling for PPV 
vertices as mentioned above: 



gKKu 



(5) 



This is the coupling constant used in the Jiilich KN mod- 
els [Tl"lll2| . In case of the D meson the coupling constant 
is given by 



9ddu 



-x9ddui 1& - 9KKp ■ (6) 



Summarizing the above results from Eqs. (|3l5l6p we see 
that 



9dd p - 9kk p 
9dDu = 9kku ■ 



(7) 



Thus, the coupling constants of the exchanged vector me- 
sons are the same for the KN and DN systems under 
assumption of SU (4) symmetry and ideal mixing. 

Since some of the couplings involving the D meson 
are known empirically, at least to some extent, we want 
to review them briefly here. The DDp coupling constant 
was determined in Refs. [19,20 based on the vector dom- 
inance model and found to be gDDp — 2.52 — 2.8. This 
value, which was subsequently adopted in several inves- 
tigations [7,21, 22J, is only marginally smaller than the 
one which follows from assuming SU(4) symmetry. The 
same is true for the DDto coupling constant, found to be 
gDDuj = —2.84 in Ref. [20 , likewise derived within the vec- 
tor dominance model. In Ref. [52] the value g^DD* — 5.56 
is cited, derived from the measured decay width of the D* 
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Table 1. SU(4) matrix representation of the pseudo-scalar (P) and vector (V) mesons. 
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Table 2. Vertex parameters used in the meson-exchange model of the DN interaction at the DDm r or DD*m r (M) and NNm r 
or NAm r (B) vertices. M r and m r refers to the mass of the exchanged particle. Note that the scalar meson exchanges (a, ao) 
are considered as an effective interaction, cf. text, and therefore we provide only the product of the coupling constants. 



Process 


Exch. part. 


M r or m r 


gM/v4n 


g B /V^ [/b/^b] 


A M 


A B 






[MeV] 






[GeV] 


[GeV] 


DN -> DN 


P 


769 


0.843 


0.917 [6.1] 


1.4 


1.6 






782.6 


0.843 


2.750 [0.0] 


1.5 


1.5 




Ac 


2285 


-2.284 


-2.284 


4.1 


4.1 






2455 


0.435 


0.435 


4.1 


4.1 




a 


600 




0.25 (1.00) 


1.7 


1.2 






980 




0.65 (2.60) 


1.5 


1.5 


DN -> D*N 


TV 


138.03 


0.843 


3.795 


1.3 


0.8 




p 


769 


0.843 


0.917 [6.1] 


1.4 


1.0 


DN -> D*A 


TV 


138.03 


0.843 


0.600 


1.2 


0.8 




P 


769 


0.843 


5.740 


1.3 


1.0 


DN -> DA 


P 


769 


0.843 


5.470 


1.3 


1.6 



meson. Here the corresponding SU(4) coupling constant 
is roughly a factor 2 smaller. 

In any case, in our model calculation we use coupling 
constants that are determined fully by SU(4) symmetry. 
The difference to those values deduced from available ex- 
perimental information is not very large and, thus, does 
not really warrant a departure from SU(4) at present. In- 
deed, there are other assumptions made in the model cal- 
culation, that could be considered to be more question- 
able, for example those about the vertex form factor. As 
mentioned above, the prescription we use relies on the 
fact that the same particles are exchanged in the KN 
and DN potentials. Possible influences from differences in 
the off-mass-shcll dependence due to the different (KN or 
Z)iV) intermediate states, appearing in higher iterations, 
are simply ignored. However, the main uncertainty in the 
meson-exchange model arises from the treatment of the 
scalar-meson sector. Here, unlike for pseudoscalar and vec- 
tor mesons, so far there is no general agreement about who 
are the actual members of the lowest lying scalar-meson 
SU(3) multiplet. (For a thorough discussion on that is- 
sue and an overview of the extensive literature we refer 
the reader to [2"51l2~4"j and references therein.) Therefore, 
it remains unclear whether and how the relations for the 



coupling constants given in Eq. (|TJ) should be applied in 
the SU(3) case [12,25], but even more so when it comes 
to SU(4). It is known for a long time that the contri- 
butions from the scalar sector play a crucial role in any 
baryon-baryon and meson-baryon interaction at interme- 
diate ranges. 

In the present paper we consider two different scenar- 
ios for the scalar mesons. First, in line with the works 
in Refs. [12 ,25], we view the contributions in the scalar 
sector as being due to correlated irir and KK exchange. 
However, in the absence of a concrete model for those 
contributions, as it was used in [T2125] , we resort here to 
a rough estimation of their strength. Based on the scale 
of the correlated tttt and KK exchange, which we iden- 
tify with the masses of the relevant propagators - p (K*) 
and/or tttt ( KK) for the ttN and KN interactions, but 
the D*(2010) and/or DD for the DN system - we ex- 
pect that its strength should be about 4 times smaller in 
the latter case. Thus, we reduce the coupling constants by 
that factor as compared to the values used in [10] . Note 
that this reduction is supported by available model cal- 
culations of the reaction pp — > DD [21)112"?] which suggest 
that the corresponding amplitude, which would form the 
main ingredient for a microscopic calculation of the cor- 
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related 7T7t exchange for DN, cf. Ref. [12], is significantly 
smaller than the one for pp — > KK, even when taking into 
account the kinematical differences. The used values are 
given in Table [2] The second scenario is an attempt to 
simulate the case that the scalar contributions are due to 
genuine scalar-meson exchange. Accordingly, we use the 
same scalar coupling constants for the DN interaction as 
in the KN model [10j . The concrete values for that sce- 
nario are listed in brackets in Table [2] 



in 




3 The DN interaction based on the 
quark-gluon exchange mechanism 

The quark interchange processes with one-gluon-exchange 
(OGE) we consider in the present paper are represented 
pictorially in Fig. In these graphs, the D mesons are 
D° = uc and D~ = dc, so that the exchanged quarks 
are always the light u and d quarks, the c antiquarks are 
not interchanged. The dominant contributions of the OGE 
interaction are the Coulomb and spin-spin parts. The in- 
teraction of the quarks i and j with constituent masses 
m.i and rrij can be written as Vy = TfT" Vij(q 7 S), where 
T a = X a /2 for a quark and T a = -(A Q ) T /2 for an anti- 
quark, and the momentum and spin dependent pieces as 
Vij (g, S) are given as 



vc(q) 

47ra, 



q 2 3mirrij 



v S s(q) Si-Sj 



(8) 



where a s is the quark-gluon coupling constant. 

As shown in more detail in Appendix (B[ the effective 
DN potential Vq n can be written as a sum of four con- 
tributions as 

4 

Vdn(p,p') = J2^ [y i (P,P') + V l (p , , P )}/2, (9) 

i=l 

where each term in the sum corresponds to a graph in 
Fig. [2 The uji's are given Tabled they come from sum- 
ming over the color-flavor-spin indices of the quarks and 
include symmetry combinatorial factors. The V^s are func- 
tions of the center-of-mass momenta p and p ' which are 
given by multidimensional overlap integrals over the in- 
ternal wave functions of the nucleons and mesons and the 
OGE potentials. When using Gaussian forms for the nu- 
cleon and meson wave functions, many of the integrals can 
be done analytically and the V^'s can then be expressed in 
terms of a single three-dimensional integral as 



Vi(p,p') = e -^P 2 ~b z p' 2 +c tP - P ' 



3.9 



[(3 + 2g)Tra 2 N 



3/2 



d 3 q 



v{q)e~ 



m +ei-q 



(10) 



.9 



where v(q) = t>c(g), or vss(q) defined in Eq. 
a 2 N j f$ 2 D , where ajv and /3d are the Gaussian widths of the 
nuclcon and the D meson. The a*, &<, • ■ • are given in terms 
of g and the quark masses - see Appendix [Bl 



(3) 



(4) 



Fig. 2. Pictorial representation of the four different quark- 
interchange processes that contribute to the DiV interaction. 
The wavy lines represent the one-gluon exchange (OGE) and 
the solid lines represent quarks. 



4 Results and discussion 

The original KN model of the Julich group includes be- 
sides the exchange of the standard mesons also an addi- 
tional phenomenological (extremely short-ranged) repul- 
sive contribution, a "oy ep ", with a mass of about 1.2 Gev 
[T2] . This contribution was introduced ad- hoc in order to 
achieve a simultaneous description of the empirical KN 
S- and P-wave phase shifts, but it is also required for a 
consistent description of the KN and KN systems |28j . 
which are interrelated via a G-Parity transformation. Ev- 
idently, due to its phenomenological nature it remains un- 
clear how that contribution should be treated when go- 
ing over to the DN system. Fortunately, a recent investi- 
gation by our group provided evidence that a significant 
part of that short-ranged repulsion required in the origi- 
nal Julich model could be due to genuine quark-gluon ex- 
change processes. Thus, in the present study we will build 
upon this insight when constructing a model of the DN 
interaction. In particular this means that also for the DN 
system we consider contributions from meson-exchange as 
well as from quark-gluon mechanisms where each of them 
is closely linked to the corresponding pieces in the KN in- 
teraction. However, in order to get a better understanding 
on what changes and what remains the same in the transi- 
tion from KN to DN we will first study those sectors sepa- 
rately. In order to facilitate an easy comparison of the KN 
and DN results we present them as a function of the cor- 
responding excess energies. (One should be aware that a 
comparison at the same laboratory momentum, say, would 
look quite different because of the large mass difference be- 
tween the kaon and the D meson.) Within the range shown 
the KN and DN reactions are predominantly elastic. The 
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Table 3. The color-spin-flavor coefficients from the spin-spin OGE interaction for the uii for the D N and D N systems for 
the individual charge states, and for the 1 = and 1=1 combined isospin states. 



Process 






U>2 






^3 




UJ4 






Si Sj 




SiSj 




SiSj 




SiSj 


pD° -*pD° 


1/3 


1/3 


1/3 


1/3 


1/3 


1/18 


1/3 


1/18 


n D~ — ► n D~ 


1/3 


1/3 


1/3 


1/3 


1/3 


1/18 


1/3 


1/18 


pD~ =>pD~ 


1/3 


1/6 


1/3 


1/6 


1/3 


1/9 


1/3 


1/9 


nD° ~*nD° 


1/3 


1/6 


1/3 


1/6 


1/3 


1/9 


1/3 


1/9 


pD~ -> nD° 


1/3 


1/6 


1/3 


1/6 


1/3 


-1/18 


1/3 


-1/18 


1 = 

















-1/6 





- 1/6 


I = 1 


- 4/9 


-1/3 


+ 4/9 


- 1/3 


+ 4/9 


- 1/18 


-4/9 


- 1/18 



first inelastic hadronic channel (pion production) opens at 
an excess energy of around 136 MeV. We want to mention 
also that, due to the much smaller mass difference between 
D*(2010) and L»(1869) versus K*(892) and X(496), the 
nominal threshold of the D*N channel occurs at a signifi- 
cantly smaller excess energy than the corresponding K*N 
channel. Indeed the former practically coincides with the 
DNir threshold. 

Based on the DN interaction potential V described 
in the two preceeding sections the corresponding reaction 
amplitude T is obtained by solving a Lippmann-Schwinger 
type scattering equation defined by the time-ordered per- 
turbation theory, 

r = v + vg r , 

from which we calculate the DN observables in the stan- 
dard way [11] . Due to the large mass of the D meson that 
enters into the free Green's function Go higher iterations 
play a somewhat less important role for DN as compared 
to the KN system. But we want to emphasize that uni- 
tarization of the reaction amplitude, which is achieved by 
solving Eq. is essential for obtaining meaningful re- 

sults because the resulting phase shifts in the S- and also 
P-waves are in the order of 20 degrees or even more in 
the energy range covered by our study. For completeness 
let us also mention that we use averaged masses for the D 
mesons, namely m£>=1866.9 MeV and m£>.=2009 MeV. 

Results based on the meson-exchange contributions are 
shown in Figs.[3](for isospin 1 = 0) and[4](for 1 = 1). Since 
some of the DN model calculations in the literature take 
into account only p exchange [7l[9] we consider its contri- 
bution first. It is obvious that the resulting cross sections 
for the DN system are much larger than those for KN. 
But one should keep in mind that this difference is pri- 
marily caused by the different kinematics (masses). The 
involved coupling constants are exactly the same for the 
KN and DN interactions under the assumption of SU(4) 
symmetry, as discussed in Sect. 2. It is worth noting that 
the cross section in the 1 = channel of the DN is partic- 
ularly large and even exceeds 100 mb near threshold, while 
for I = 1 is is less then 5 mb over the whole considered 
energy range. 

The picture changes drastically once uj exchange is 
added. Specifically, in the 1 = channel its contribution 
interferes distructively with the (attractive) p exchange 
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Fig. 3. KN and DN cross sections in the isospin channel 
1=0 including consecutively p (dashed curve), uj (dash-dotted), 
scalar mesons and baryon-exchange diagrams (dotted), and 
box diagrams (solid). 
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Fig. 4. KN and DN cross sections in the isospin channel 1=1. 
Same description of curves as in Fig. [3] 



and leads to a strong reduction of the predicted cross sec- 
tion. On the other hand, for 7=1 both contributions are 
repulsive and add up so that now the cross sections in 
both isospin channels are of comparable magnitude. In- 
deed after inclusion of the w exchange the KN results 
show already the typical features of the full model [12] 
but also of the experimental information [29j . namely an 
almost constant cross section for 1 = 1 and a cross section 
for 1 = that is practically zero at threshold and then 
increases with energy. The predictions for the DN system 
exhibit very similar features. 

The addition of the scalar contributions and of baryon 
(yt c (2285), I7 C (2455)) exchange influences the results for 
DN very little and therefore we don't show them sepa- 
rately. This is not too surprising in view of our assump- 
tion about the origin of the scalar sector (we will come 
back to this issue later) and of the large mass of the ex- 
changed baryons. On the other hand, the box diagrams 
yield a sizeable contribution, in particular in the 1 = 
channel of the DN system. 



J2 

E 
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Fig. 5. KN and DN cross sections from OGE. The dashed- 
dotted (dashed) curve is the result for the spin-spin (SS) part 
in the / = 1 (I = 0) channel. The solid curve is the result for 
the 1 = 1 channel after adding the Coulomb part. The dotted 
curve is the corresponding cross section as obtained in Born 
approximation. For the 1 = channel the Coulomb component 
is zero so that the full OGE result coincides with the curve for 
the SS part. 



Results for the KN and DN interactions based on 
the OGE in the quark model are presented in Fig. [Sj We 
use standard quark model parameters |30j . For the light 
quark masses we take m u = md = 330 MeV and for the 
strange and charm quark masses we use m s = 550 MeV 
and m c = 1600 MeV, and the quark-gluon coupling is 
taken to be a s = 0.6. The size parameters of the nu- 
cleon and the kaon wave functions are ajv = 400 MeV 
and (3d = 350 MeV. For the D meson wave function we 
use the value of Ref. [31], namely (3 D = 383.5 MeV. The 
results demonstrate that, as in the case of KN, the spin- 
spin component of the OGE is much more important than 
the Coulomb component. (Note that the contribution of 
the Coulomb component is zero in the 1 = channel.) 
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We have performed exploratory calculations utilizing a 
larger value for (3d for the D meson, (3d — 440 MeV, 
as given by a recent calculation [52] . The results for the 
combined Coulomb and spin-spin OGE do not change ap- 
preciably, although the Coulomb part is a little smaller 
in this case as compared to the corresponding value with 
(3 D = 383.5 MeV. 

The cross sections predicted for the DN system are 
roughly a factor 1.5 larger than those for KN. Note that, 
unlike for the meson-exchange part, here the parameters 
entering the potential differ, reflecting the different quark 
masses and sizes of the K and D mesons. As in the KN 
case, graphs (2) and (4) in Fig. [5J are suppressed as com- 
pared to graphs (1) and (3), because of the large quark 
mass in the denominator of Eq. ([8]). This suppression is 
even stronger in the DN case. Note that graphs (1) and 
(3) actually become somewhat larger as compared to the 
kaon case because of the larger size of the D meson. 

It is worth pointing out there are sizable effects due to 
the iteration of the interaction in the Lippmann-Schwinger 
equation. In order to demonstrate this we include also re- 
sults obtained in Born approximation for the 1=1 chan- 
nel. These are shown by the dotted lines in Fig. [51 Obvi- 
ously, close to threshold the Born result for the DN cross 
section is of the order of 30 mb, while the corresponding 
unitarized result is of the order of 10 mb, i.e. there is more 
than 50 % difference. 

Our full predictions for the DN cross sections, com- 
bining now the mesonic part with the contributions from 
quark-gluon processes, are presented in Fig. [6l The solid 
lines are results for the scenario where the scalar-meson 
contributions are viewed as being due to correlated tttt — 
KK exchange, in line with the philosophy of the original 
Jiilich KN model [12]. In this case the cross section for 
J = 1 is in the order of 20 mb, i.e. roughly twice as large 
as observed for the KN system [29]. For the 1 = chan- 
nel we predict a cross section that is practically zero at 
the threshold but increases to about 25 mb at the excess 
energy 150 MeV. Also here the result is roughly twice as 
large as the cross section for KN at the corresponding 
excess energy. 

The dashed curves show results obtained in the sce- 
nario that attempts to simulate the case that the scalar 
contributions are due to genuine scalar-meson exchange. 
In order to get a rough estimate for that scenario we as- 
sume here, for simplicity reasons, that the couplings for 
DN are the same as for KN for the two scalar mesons 
in question. A strict evaluation within the SU(4) scheme 
would involve several unknown quantities such as the sin- 
glet couplings and the mixing angles and is not feasible. 
In any case, our prescription is only meant to illustrate 
the ambiguity resulting from the unclear situation in the 
scalar sector. Evidently, the cross section for J = is 
rather insensitive to the treatment of the scalar mesons, 
at least within the scenarios considered here. Indeed, due 
to their isospin structure, the contributions of the a and 
the do (980) mesons tend to cancel in this channel. On the 
other hand, the predicted 1=1 cross section reduces by 
50 % for the scenario based on larger scalar-meson cou- 
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Fig. 6. DN cross sections in the isospin channels 1=0,1. The 
solid curves are the results of the full model, i.e. including 
meson-exchange and OGE, and viewing the scalar contribu- 
tions to be due to correlated rnv exchange. The dashed curves 
show results obtained in the scenario that attempts to simu- 
late the case that the scalar contributions are due to genuine 
scalar-meson exchange, cf. text. 



pling constants. This variation may be considered as a 
measure for the uncertainty in our model prediction for 
DN, despite of constructing the interaction in close anal- 
ogy to KN and invoking strict SU(4) symmetry. 

Finally, for completeness let us also present the S- 
wave DN phase shifts. Corresponding results are shown 
in Fig. [3 for the quark-gluon interaction alone (upper 
panel) and for the full model (lower panel). Obviously, 
in general all interactions are repulsive, as reflected in the 
negative sign of the phase shifts. But the full model is 
weakly attractive for energies near the threshold in the 
1 = channel. The corresponding scattering lengths are 
a 7=0 = — 0.13fm, a I=1 = — 0.29fm, for the quark-gluon 
interaction and a I=0 = 0.07 fm, a /=1 = — 0.45fm, for the 
full model. Interestingly, the former results are pretty close 
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Fig. 7. DN 5-wave phase shifts in the isospin channels 1=0,1. 
The upper panel shows the results for the quark model based 
on OGE while the lower pannel are the results of the full model, 
i.e. including meson-exchange and OGE, and viewing the scalar 
contributions to be due to correlated mr exchange. 

to the values found by Lutz and Korpa for their DN inter- 
action [16] while the latter are qualitatively very similar 
to the results obtained for the KN interaction [12] . 

5 Summary 

In this paper we presented predictions for the low-energy 
DN cross section based on a model which was developed 
in close analogy to the meson-exchange KN interaction 
of the Jiilich group [TllfT2] . utilizing SU(4) symmetry con- 
straints. The main ingredients of the interaction are pro- 
vided by vector meson (p, u>) exchange but higher-order 
box diagrams involving D*N, DA, and D* A intermedi- 
ate states, are taken into account too. Furthermore, in the 
spirit of a recent study of the KN system by us [10] , the 
short range part is again assumed to receive additional 
contributions from genuine quark-gluon processes. 



The cross sections for DN predicted for excess ener- 
gies up to 150 MeV are of the same order of magnitude 
as those for KN but with average values of around 20 mb 
roughly a factor two larger than for the latter system. 
There is an uncertainty in our prediction for the 1 = 1 
channel which is caused by the unknown SU(4) structure 
of the scalar-meson sector. Assuming that the contribu- 
tions in the scalar sector are due to correlated 7T7t ex- 
change, in line with the Jiilich KN model, we find that 
the scalar contributions influence the DN cross sections 
only marginally. But a scenario where the effect of the 
exchange of genuine scalar mesons is simulated by assum- 
ing their coupling strengths to be the same as in the KN 
model leads to a 50 % reduction of the 1=1 cross section. 

Anyway, the most interesting finding of our study is 
certainly the important role played by the w-exchange 
contribution. Its interference pattern with the p-exchange, 
which is basically fixed by the assumed SU(4) symmetry, 
clearly determines the qualitative features of the DN in- 
teraction - similiar to what happens also for the KN sys- 
tem^ As a consequence also the cross sections predicted 
for DN show qualitatively very similar features to those 
known of KN scattering. On the other hand, predictions 
for DN where only p-exchange was taken into account dif- 
fer drastically and, in our opinon, should be regarded with 
caution in view of the results presented in this paper. 

Acknowledgements 

This work was financially supported by the Deutsche For- 
schungsgemeinschaft (Project no. 444 BRA-113/14) and 
the Brazilian agencies CAPES, CNPq and FAPESP. This 
work was supported in part by the EU I3HP "Study of 
Strongly Interacting Matter" under contract number RII3- 
CT-2004-506078, by the DFG through funds provided to 
the SFB/TR 16 "Subnuclear Structure of Matter", by the 
EU Contract No. MRTN-CT-2006-035482, "FLAVIAnet" 
and by BMBF (grant 06BN411). A.S. acknowledges sup- 
port by the JLab grant SURA-06-C0452 and the COSY 
FFE grant No. 41760632 (COSY-085). 



10 



J. Haidenbauer et al.: DN interaction 



A The interaction Lagrangians 

Here we list the specific interaction Lagrangians which are 
used to derived the meson-exchange DN interaction. The 
baryon-baryon-meson couplings are given by 



£nns 
Cnnp 
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9NNP&N {x)rf#N {x)$P (X) 
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■(&y(x)'y b 'fW N (x) 



d v $ v (x)) 



(11) 



Here <£at, ^/V; an d are the nucleon, A, and hyperon 
field operators and <Ps, ^p> an d are tne ne ld operators 
for scalar, pseudoscalar and vector mesons, respectively. 
The employed three-meson couplings are 



£-pps 
Cppv 

Cvvp 



gpp S m P <P P (x)<P P (x)<P s (x) 
g PP v$p{x)d^p(x)®v{x) 



gwp 



ie^ TS d^» v (x)d T $ s v {x)$p{x) , 



where t^ VT s is the antisymmetric tensor with e 



0123 



(12) 
= 1. 



Note that here only the space-spin part is given. The addi- 
tional SU(4) flavour structure that leads to the character- 
istic relations between the coupling constants is discussed 
in Sect. II. Details on the derivation of the meson-baryon 
interaction potential from those Lagrangians can be found 
in Refs. QJIE2]. 



B The quark-model meson-baryon interaction 

In this Appendix we outline the derivation of the effective 
meson-baryon interaction in the quark model. As already 
mentioned, the effective interaction is given by the quark- 
Born diagrams depicted in Fig. [2j Their expressions can be 
obtained in the Born-order quark interchange model |33] , 
or using the methods of the resonating group, or the Fock- 
Tani representation [34]. Given the interactions and the 
bound state amplitudes of the single hadrons, the expres- 
sion of the quark Born diagrams for the effective meson- 
baryon interaction a + [3 —> "f + 5 is given by 



V 



MB 



V qq {pv-ap)^ v ^ 
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3 </>•" /' ij/' 1 ' " 



6C 



V qq {pv;ap)^"^ 



W2W 
CT/i3 



Here, the $ and \P are Fock-space amplitudes of the one- 
meson and one-baryon states, which in a second quanti- 
zation notation are given as 



M a ) = $Z»qUi\0), \B a ) = 



1 



)B/*l/*I/i3 „t fl t J 

a y^i Hfj, 3 



|0>, 
(14) 

where a indicates all quantum numbers necessary to spec- 
ify the hadronic state, like cm. momentum, spin and fla- 
vor, and p, v, ■ ■ ■ indicate all quantum numbers of the 
quarks like momentum, color, spin and flavor - a sum 
or integral over repeated indices is implied. q> , q\ q\ and 
q are quark and antiquark creation and anhihilation op- 
erators that satisfy the usual canonical anticommutation 
relations. In addition, V qq , V qq and V qq are the microscopic 
quark and antiquark interactions, which in the same sec- 
ond quantization notation are defined through 



1, 



y = ^Vqqipv'i P^hl ll <1<j 1 P + 7;V<iq(V>v, po)ql ql q„ q P 



t /7t 



V qq (pu; pa)q [ql q f a q 



(15) 



The expression for Vmb in Eq. (H2J) involves a 6-dimen- 
sional integral that cannot be integrated analytically for 
general forms of the amplitudes and W. However, when 
using Gaussian forms for the meson and baryon ampli- 
tudes, many of the integrals can be done analytically and 
the resulting expression for each of the diagrams of Fig. [5] 
is of the form given in Eq. (fT0|) . Specifically, we use for the 
amplitudes <P and )P in momentum space Gaussian forms 
with width parameters (3d a nd Q-n as 



<Mfci,fc 2 ) = <5 (3) [P-kx-k, 
where 



7T/3* 



3/2 



2(M g fci - M q k 2 ) 
M q + M q 
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(17) 



and 



V P (ki,k 2 ,k 3 )= sw(p 
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V qq (pv;ap) 



k p )v(k^ 



(18) 



(19) 



and equivalently for V qq and V qq , after integrating over the 
quark momenta as indicated in Eq. (|13p one obtains the 
expression given in Eq. (1101 ), where the a^, hi, ■ ■ ■ , can be 



written as a ratio di = n(ai) / ' d(a,i)a 2 Nl bi — n(bi) / d(bi)c 
(13) etc. The corresponding expressions are given as follows: 



■N' 
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Graph (1): 



n(oi) = 3.g 2 
d(ai) = 6(3 



- 3 (1 + p) 2 + .g (7 
2g) (l + p) 2 
n(ci) = .g 2 + (l + p) 2 -2 5 ' 



10 p 2 



d{c x ) = (3 + 2.g) (1 + p) 2 
n(di) = 3g, d(ai) = (3 + 
n(ei) = -g (l+4p) (p + p 
d( ei ) = (3 + 2 ff ) (1+p), 



-1 + ^ 



2g) 



Graph (2): 



Graph (3): 

n(oa) - 3.g 2 + 3 (1 + p) 2 + .g (7 + 8p + 10p 2 ) 

d(a 3 ) = 6 (3 + 2 5 ) (1 + p) 2 

"(63) = 1(03), d(b 3 ) = d(a 3 ) 

n(c 3 ) = n(ci), d(c 3 ) = d(ci) 

n(d 3 ) = 6 + 7 5 , d(d 3 ) =4 (3 + 2.g) 

n(e 3 ) = -3 [1 + 5 + (1 + 2.g) p] p + [3 + .g + (3 

d(e 3 ) = 2 (3 + 2 5 ) (1 + p), 



25) P]p' 
(22) 



Graph (4): 



In these, g = a%/f3p and p = M u /M c 



References 

1. T. Matsui and H. Satz, Phys. Lett. B 178, 416 (1986). 

2. R. L. Thews, M. S chroedter and J. Rafelsk i, Phys. Rev. C 
63, 054905 (2001) arXiv:hep-ph/0007323j . 

3. H. Satz, J. Ph ys. G 32, R25 (2006) 
[arXiv:hep-ph/0512217] . 

al., Technical Progress Re- 
Strong Interaction Stud- 
February 2005, 



M. Kotulla et 
port for PANDA, 
ies with Antiprotons, 



10. 



11. 



12. 



13. 



n(a 2 ) 


= n(oi), d(a 2 ) = d(ai) 






n(b 2 ) 


= n(6i), d(6 2 ) = d(h) 




14 


n(c 2 ) 


= n(ci), d(c 2 ) = d(ci) 




15 


n(d 2 ) 


= ff (3 + 5 ), d(d 2 ) = 2(3- 


f 2 5 ) 




n(e 2 ) 


= g [(2 + .g + 2p)p-(l + .g 


-2p)p'] 


16 


d(e 2 ) 


= (3 + 2 5 ) (1 + p), 




( 21 ) 17 



71(04) 


= n(ai), 


71(64) = 


n(6 x ) 


25 


d(d4) 


= d(ai), 


d(6 4 ) = 


d(h) 


n(c 4 ) 


= n(a), 


d(c 4 ) = 


d(ci) 


26 


n(d4) 


= "2 + 9, 


d(d 4 ) = 


4 


27 


n(e 4 ) 


= -(1 + 


.9 + P) (P - 


-pO 


d(e 4 ) 


= 2 (1 + 


P). 




(23) 28 



W. Cassing, Y. S. Golubeva and L. A. Kondratyuk, Eur. 
Phys. J. A 7, 279 (2000) [arXiv:nucl-th/9911026] . 
T. Barnes, |arXiv:nucl-tn/030603ir ~ 

Z. w. Li n, C. M. Ko and B. Zhan g, Phys. Rev. C 61, 024904 

(2000) [arXivmucl th/9905003] . 

A. Sibirtsev, K. Tsushima and A. W. Thomas, Eur. Phys. 
J. A 6, 351 (1999) [arXiv:nucl-th/9904016| . 
A. Sibirtsev, Nucl. Phys. A 680, 274c (2001). 

D. Hadjimichef, J. Haidenbauer and G. Krein, Phys. Rev. 
C 66, 055214 (2002) [arXivinucl-th /0209026] . 
R. Biittgen, K. Holinde, A. Miiller-Groeling, J. Speth, and 

(20) p. Wyborny, Nucl. Phys. A506, 586 (1990). 

M. Hoffmann, J.W. Durso, K. Holinde, B.C. Pearce, and 
J. Speth, Nucl. Phys. A593, 341 (1995). 
J. Hai denbauer and G. Krein, Phys. Rev. C 68, 052201 
(2003) [arXiv:hep-ph /0309243| . 

J. Hof mann and M. F. M. Lu tz, Nucl. Phys. A 763, 90 

(2005) [arXiv:hep-ph/0507071| . 

L. Tolos, J. Schaffner-Bielich and A. Mishra, Phys. Rev. C 
70, 025203 (2004) [arXiv:nucl-th/0404064| . 
M. F. M. Lutz and C. L. Kor pa, Phys. Lett. B 633, 43 

(2006) [arXiv:nucl-th/0510006l . 

T. Miz utani and A. Ramos, Phys. Rev. C 74, 065201 
(2006) [arXiv:hep-ph/0607257l . 

18. R. Machleidt, K. Holinde, Ch. Elster, Phys. Rep. 149, 1 
(1987). 

19. S.G. Matinyan and B. Miiller, Phys. Rev. C 58, 2994 
(1998). Note that the actual coupling constant given in 
that reference is a factor of 2 larger due to the difference 
of the definitions. 

20. Z.-w. Lin and C. M. Ko , Phys. Rev. C 62, 034903 (2000) 
[arXiv:nucl-th/9912046] . 

21. Z.-w. L in, T. G. Pi and C. M. K o, Nucl. Phys. A 689, 965 

(2001) [arXiv:nucl-th/0006086l . 

W. Liu and C. M. Ko , Phys. Lett. B 533, 259 (2002) 
|arXiv:nucl-th/0201074]. 

E. Klempt, |arXiv:hep-ph/0404270l 

Y. Kalashnikova, A. E. Kudryavtsev, A. V. Nefediev, 
J. Haid enbauer, and C. Hanhar t, Phys. Rev. C 73, 045203 
(2006) [arXiv:nucl-th/0512028l . 

J. Haidenbaue r and U.-G. MeiBner, Phys. Rev. C 72, 
044005 (2005) [arXiv:nucl-th/0506019j . 
P. Kroll, B. Quadder and W. Schweiger, Nucl. Phys. B 
316, 373 (1989). 

A. B. Kaidalov and P. E. Volkovitsky, Z. Phys. C 63, 517 
(1994). 

A. Miiller-Groeling, K. Holinde, and J. Speth, Nucl. Phys. 
A513, 557 (1990). 

A. Sibirtsev, J. Haidenbauer, S. Krewald and 
U.-G. Meifiner, J. Phys. G 32, R395 (2006) 
[arXiv:nucl-th/0608028] . 

T. Barnes and E. S. Swanson, Phys. Rev. C 49, 1166 
(1994). 

K. Martins, D. Blaschke and E. Quack, Phys. Rev. C 51, 
2723 (1995) [arXiv:hep-ph/9411302j . 

J. P. Hilbert, N. Black, T. Barnes and E. S. Swanson, 
|arXiv:nucl-th /0701087| 

T. Barnes, E. S. Swanson and J. D. Weinstein, Phys. Rev. 
D 46, 4868 (1992) [arXiv:hep-ph/9207251] . 
D. Hadjimichef, G. Krein, S. Szpigel and J. S. da Veiga, 



22. 

23. 
24. 



29. 



30. 



31. 



32. 



33. 



34. 



Annals Phys. 268, 105 (1998) arXiv:hep-ph/9805459 



http:/ /www-panda.gsi.de/db/papersDB/PC19-050217_panda_tpr.pdf 



